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ABSTRACT 

Understanding gas flows into and out of the most massive dark matter structures in 
our Universe, galaxy clusters, is fundamental to understanding their evolution. Gas in 
clusters is well studied in the hot (> 10 6 K) and cold (< 10 4 K) regimes, but the warm 
gas component (10 4 - 10 6 K) is poorly constrained. It is challenging to observe directly, 
but can be probed through Lya absorption studies. We produce the first systematic 
study of the warm gas content of galaxy clusters through synthetic Lya absorption 
studies using cosmological simulations of two galaxy clusters produced with Enzo. We 
explore the spatial and kinematic properties of our cluster absorbers, and show that the 
majority of the identified absorbers are due to fast moving gas associated with cluster 
infall from IGM filaments. Towards the center of the clusters, however, the warm 
IGM filaments are no longer dominant and the absorbers tend to have higher column 
densities and metallicities, representing stripped galaxy material. We predict that the 
absorber velocity distribution should generally be bi-modal and discuss the effects 
of cluster size, mass, and morphology on the properties of the identified absorbers, 
and the overall cluster warm gas content. We find tentative evidence for a change in 
the well known increasing Ahi with decreasing impact parameter for the most massive 
dark matter halos. Our results are compared directly to observations of Lya absorbers 
in the Virgo cluster, and provide predictions for future Lya absorption studies. 
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1 INTRODUCTION 

Galaxy clusters are the largest virialized structures in our 
universe. Clusters host dark matter potential wells around 
10 14 -10 15 Mq that are able to heat the infalling gas to 
temperatures of 10'-10 S K. This hot intracluster medium 
(ICM), constitutes the most significant baryonic component 
of galaxy clusters and is well observed in thermal X-ray emis¬ 


sion (Voit 2005 and references therein). The flow of gas 


into and out of galaxies within clusters, and into the cluster 
itself, is strongly affected by the cluster environment. Our 
understanding of the competing processes involved is still 
forming, from galaxy level processes such as galactic out¬ 
flows and inflows to cluster specific processes such as ram- 
pressure stripping, galaxy mergers, galaxy harassment, and 


starvation (e.g., Cayatte et al. 1990 Vollmer et al. 2001 


Lewis et al., 2002; Roediger & Hensler[ |2005| 

Chung et al. 

2007 

Tonnesen, Bryan & van Gorkom, 2007 

Yoon et al. 


2012). In addition, our understanding of gas flows onto clus¬ 


ters from IGM filaments is incomplete. Much of this gas is 
expected to exist in a warm, diffuse phase at temperatures 


e-mail: emerick@astro.columbia.edu 


near T ~ 10 4 - 10° K (e.g., Burns, Skillman & O’Shea 2010 


Werner et al. 2008 Kravtsov & Borgani 2012). Although 


much of the baryonic mass in the Universe exists in this 
phase, it remains relatively unexplored in a cluster environ¬ 
ment. This is partly because this regime is too diffuse to 
observe directly via 21cm emission, and too cool to be seen 
in X-ray. However, Lya absorption along lines of sight to¬ 
wards background QSO’s presents a valuable means to probe 


the warm gas (e.g., Tripp et al. 2005). 

Lya absorption along lines of sight towards background 
QSO’s has long been a useful probe of the intergalactic 
medium (IGM) at high redshifts (z ~ 2 - 6) (see Rauch 
(1990 and references therein). A majority of the baryons 
in the Universe reside in a diffuse state that comprises the 
IGM. Lya absorption can account for nearly all of the high 
redshift baryons in the Universe ( |Weinberg et al (j |1997[ ). 
Our knowledge of the Lya forest in the nearby Universe is 
comparatively lacking, although attempts have been made 
recently towards understanding this redshift regime both ob- 


servationally (e.g., Danforth et al.| 20141 and with simula¬ 
tions (e.g., Dave et al. 2010 Kollmeier et al.| 2014 1 . Lya ab¬ 
sorption at z < 2 is possible with UV spectrographs aboard 
HST — most recently with the Cosmic Origins Spectrograph 
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(COS). Although Lya has often been used to study the 


IGM and warm gas associated with galaxies (e.g., Tumlin- 


son et al.l 2013), there has been little examination of warm, 


Lya absorbers in the cluster environment. To date |Yoon| 
et al. (2012) (hereafter Y12) has been the only systematic 
study of warm gas via Lya absorption in galaxy clusters. 
Y12 examined 43 Lya absorbers along 23 sightlines towards 
the Virgo cluster. They found that the warm gas in Virgo 
is preferentially located in the cluster outskirts and associ¬ 
ated with substructures. They conclude that the observed 
gas primarily represents gas flow into the Virgo cluster. To 
extend the observations of Y12 to morphologically different 
galaxy clusters, there are ongoing observations of Lya ab¬ 
sorption along 9 sightlines towards the Coma cluster (Yoon 
et al., in preparation) and planned observations of single 
sightlines through 11 clusters (PI: Tejos, HST PID:13833). 

The Virgo and Coma clusters provide a unique labora¬ 
tory to probe the warm gas content in two distinct clusters. 
The Virgo cluster is the closest cluster to our Milky Way, at 


16.5 Mpc (Mei et al. 2007). It is an irregular, cool-core clus¬ 


ter comprised of two distinct primary substructures, with a 


few associated minor substructures (Binggeli, Tammann & 
Sandagc, 1987a). The Virgo cluster has been well observed at 


all wavelengths (e.g., Sandage, Binggeli & Tammann, 1985 


Boselli, Gavazzi fe Sanvito| 2003 Urban et al. 201l| . This 

cluster has a virial mass and radius of 1.5 - 6.0 x 10 14 Mq 
and 1.6 Mpc respectively (see Section [2.1| for our dehnition 
of these quantities and Table [I] for references). The Coma 
cluster, located at z = 0.023, is larger and more massive 
than the Virgo cluster, with a virial mass and radius of 1.3 
- 1.4 x 10 15 Mq and 2.9 Mpc. The Coma cluster exhibits 
clear infalling substructure, as indicated by the dynamics of 


the cluster galaxies (Adami et al. 


of a radio relic (Brown & Rudnick 


20051, and the presence 


2011). The comparison 


between these two clusters allows for a better determination 
of how cluster mass and morphology affects its warm gas 
content. 

In this paper, we use the adaptive-mesh refinement code 


Enzo (Bryan et al. 20141 to simulate a Virgo like and a 


Coma like galaxy cluster, with corresponding masses and 
virial radii (R v ir) of 3.3 x 10 14 Mg and 14.5 x 10 14 Mg, 
and 1.85 Mpc and 3.0 Mpc. The physics incorporated in 
these cosmological simulations include the metagalactic ion¬ 
izing background, primordial and metal radiative cooling, a 
primordial chemistry network, and star formation and feed¬ 


back. Using the yt analysis toolkit (Turk et al. 2011), we 


observe Lya absorption via 2500 sight lines in each cluster. 
Generating synthetic spectra accounts for the observability 
of the HI gas by including thermal Doppler broadening ef¬ 
fects, gas line of sight velocities, and Gaussian random noise. 
We have developed a method to automatically identify ab¬ 
sorption features and fit Voigt-profiles adapted from previ¬ 


ous works (Egan et al. 2014 Dave et al.l 1997 Lanzetta, 


Turnshek & Wolfe 19871; this method performs well for 


weak to saturated lines, and even for the occasional some¬ 
what blended feature found at z = 0. We analyze our results 
with the goal of comparing to observables in the Virgo and 
Coma clusters. In addition, we develop an understanding of 
the properties and environment of the gas seen in absorption 
that can only be attained through simulation. 

This paper is structured as follows. In Section [2] we dis¬ 
cuss Enzo and our simulation setup, the cluster properties, 


and our synthetic spectra. We present our results in terms of 
understanding the warm gas distribution in galaxy clusters 
in Section[3] In Section[4]we compare our simulated clusters 
directly to the Virgo cluster observations in Y12 and pro¬ 
vide predictions for future observations. Finally, we discuss 
the results in Section [5] 


2 COMPUTATIONAL METHODS 

We outline the computational methods used for our study 
here, beginning with an overview of the simulation code we 
used to produce our clusters, Enzo, in Sec. 13 along with 
the included physics. We discuss the properties of our sim¬ 
ulated galaxy clusters in Sec. |2.1.1[ and compare each to 
the Virgo and Coma clusters. Our synthetic observations 
are discussed in Sec. |2.2[ from the generation of the spec¬ 
tra in Sec. |2.2.T| to the automatic fitting procedure used to 
identify absorbers in Sec. |2.2.2| 


2.1 Enzo Simulation 


We used the adaptive mesh refinement hydrodynamics + 


N-body code Enzo (Bryan et al. 2014 1 to produce cosmo¬ 


logical simulations of two separate galaxy clusters in order 
to compare the effect of cluster mass and morphology on the 
observed Lya absorption. The simulated clusters were first 
found in a low resolution, fixed grid simulation to approx¬ 
imately match the mass and size of the Virgo and Coma 
clusters. Once selected, we performed a zoom-in simulation 
on each galaxy cluster in order to achieve high spatial and 
mass resolution in the region of interest (the galaxy clus¬ 
ters). We refer to our simulated clusters as the Virgo-like 
and Coma-like clusters throughout this paper. The proper¬ 
ties of these two simulations are outlined below. 

Our Virgo-like and Coma-like galaxy clusters are simu¬ 
lated in periodic boxes 80 h ~ 1 Mpc and 160 h -1 Mpc on a 
side, respectively. We use cosmological parameters adopted 
from the WMAP5 ACDM results ( Hinshaw et al.| |2009 |. 
Specifically, fi m = 0.314, Ga = 0.686, h = Ho/100 km s _1 
Mpc -1 = 0.671, fib = 0.049, n = 0.962, and as = 0.9. The 
initial conditions are generated at z = 50 with the |Eisen-| 
stein & Hu (1999) power spectrum. The zoom-in volumes 


for each cluster were constructed to have similar spatial and 
mass resolutions, with a maximum spatial resolution for each 
cluster at z = 0 of 7.3 kpc and a dark matter particle mass 
resolution of 4.2xl0 8 Mg. 

The galaxy cluster halos were identified using a publicly 
available version of the HOP halo finding algorith m (|Eisen-| 
stein & Hut 1998 ) implemented in yt (Turk et al. 20111. In 


this paper, we use R VII to define the virial radius of a galaxy 
cluster, such that the mean density inside this radius is 100 
times the critical density of the universe at that redshift, 


consistent with the definitions in Bryan & Norman (19981. 


The HI fractions used to generate the Lya absorp¬ 
tion are obtained directly from the simulation, as we 
self-consistently solve a chemical reaction network of nine 
species: e _ , H, H + , He, He + , He ++ , H“, H 2 , and Hj (Abel 


et al., 19971. We use primordial radiative cooling based upon 


the non-equilibrium chemistry, and include metal line cool¬ 
ing as implemented in |Smith, Sigurdsson fe Abel| ( |2008| ); 
Smith et al. (20111. In solving the chemical rate equations, 
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Virgo 


ref. 


Virgo-like Coma 


ref. 


Rvir(Mpc) 

Virial Mass (10 14 Mg) 
T (keV) 

0 V (km s _1 ) 


1.6 
2.2 
~ 2.5 
544 


Mamon et al. 


iVlamon et al. 


Shibata et ai. 


2004) 


(2004) 


2001 


IVLei et ai. 12007 \ 


1.85 

3.3 

2.0 

517 


2.9 

14 

~ 8.5 
1082 


Lokas & Mamon 

(2003 

Lokas & Mamon 

(2003) 

Arnaud et al. i 

2001 

Colless & Dunn 

199BI 


Coma-like 


3.0 

14.5 

6.4 

972 


Table 1. The virial radius, virial mass, ICM gas temperature, and velocity dispersion of the Virgo and Coma galaxy clusters in comparison 
with the simulated Virgo-like and Coma-like clusters. 


we adopt the UV background from jHaardt fe Madau||2012[ ). 
We discuss the influence of our choice of ionizing background 
on our results in Sec. 15.3.11 


Our simulations include a method of star formation 
and thermal feedback which contains a simple prescription 
for purely thermal stellar feedback. The method and tests 
are described in detail elsewhere (Bryan et al.| |2014| |Smith 


et al., 2011), but we briefly summarize it here. A star forma¬ 


tion rate is computed in each cell as psfr = £SFp/idyn where 
tdyn is the local dynamical time and esf = 0.02 is the star- 
formation efficiency per free-fall time. Once created, star 
particles inject gas, thermal energy and metals into their 
local cell over a time-scale of a few dynamical times. We 
assume 10 51 erg per 55 solar masses of stars created, and a 
yield of 2%. 


Although we expect to be able to sufficiently capture the 
relevant physics for studying the warm and hot gas content 
of our clusters, we note that our simulations do not contain 
detailed physical processes that may be involved in the pro¬ 
duction or destruction of cold, dense gas in galaxies. In par¬ 
ticular, we do not expect to be able to accurately reproduce 
observed stellar and gas properties of galaxies in our simula¬ 
tions. This is because: (i) our maximum resolution is limited 
to 7.3 kpc; (ii) the thermal feedback model we adopt is rec¬ 
ognized to rapidly radiate the injected supernovae feedback 
and so fails to regulate star formation or generate galactic 
winds, and (iii) we do not include AGN feedback. The pos¬ 
sible implications of some of these processes to our work will 
be discussed in greater detail in Sec. |5.3| but to look ahead, 
this will primarily affect our results for higher column den¬ 
sity absorbers. 


We also note that we are unable use the simulations 
used in Y12 for comparisons as the simulation data used in 
that paper was corrupted in storage and is now lost. How¬ 
ever, a recent reappraisal of the existing data from Y12 
showed that the original analysis was incorrectly carried 
out on a 2 = 0.56 data set, rather than z = 0. This, 
combined with the use of a post-processing step that as¬ 
sumed a z = 0 ionizing background, would have resulted 
in an over prediction of the column densities by a factor of 
approximately (1 + 0.56) 51 ~ 10. This adjustment would 
bring the simulation results of that paper in approximate 
agreement with those found here. In addition to correcting 
this error, our current simulations are an improvement over 
that previous work since we track non-equilibrium effects on 
the HI fraction directly, while Y12 assumed ionization equi¬ 
librium and calculated the HI fraction via post processing 
with CLOUDY (Ferland et al. 1998| 20131. In addition, we 
implement the most up to date ionizing background from 
Haardt & Madau| |2012| ), which was not available at the 
time the Y12 simulations were made. 


2.1.1 Simulated Cluster Properties: Comparison to Virgo 
and Coma 


Table [l] provides a comparison of cluster properties between 
our simulated clusters and the Virgo and Coma clusters. Our 
Virgo-like cluster is slightly more massive than the Virgo 
cluster mass adopted in Mamon et al. (20041, but lies in the 


middle of the range 1.5 - 6.0x10 Mg given in Bohringer 


[et al. ( [ 1994[ >. The virial radius of the simulated cluster is 
slightly larger than that of the Virgo cluster. We therefore 
consider our Virgo-like cluster to be a somewhat larger and 
more massive version of the Virgo cluster. However, the most 
important similarities are morphological. The Virgo cluster 
is highly irregular, with several dynamically significant and 
currently merging substructures that have been well stud¬ 
ied (Binggeli, Tammann & Sandage 1987b). As shown in 


7b|. 

' 0 ° 


the density projections (top left) in Figure |H our Virgo-like 
galaxy cluster has multiple obvious and infalling substruc¬ 
tures. This cluster is then dynamically young. 

Our Coma-like galaxy cluster also has similar mass to 
(but slightly larger than) the Coma cluster. As shown in 
Figure [I] (top right), its nearly spherical shape offers a con¬ 
trast to the irregular Virgo-fike cluster. Although our Coma¬ 
like cluster is comparable in mass to Coma, it does have 
some morphological differences. Our Coma-like cluster has 
no obvious infalling substructure, yet the Coma cluster has 
a well observed infalling galaxy group, and possibly other 
substructures (e.g., Briel^Ienry_&Boehringer|^ 1992^ 


Col- 


less fe Dunn| |1996b| |Adami et al.[ |2005[ |Neumann et al. 


2003). 


The middle two panels for Figure [l] give the density 
weighted temperature projections for the two clusters. It 
is apparent that the ICM in the more massive Coma-like 
cluster is hotter, with temperatures well above 10 ' K out 
to Rvir- The temperature distribution in each cluster tends 
to follow the density projection, as the Coma-like cluster is 
more spherical than the Virgo-like cluster, and the substruc¬ 
tures seen in density in the Virgo-like cluster are apparent 
as temperature enhancements. The density weighted projec¬ 
tions allow for easy visual identification of galaxies in each 
cluster. The filaments feeding the Coma-fike cluster are ob¬ 
vious also in temperature projection, but are obscured in 
the still merging Virgo-like cluster. The lower two panels 
show the integrated HI column density projection for each 
galaxy cluster. As expected, the largest concentrations of HI 
in each cluster correspond directly to a galaxy in the given 
cluster, or tails of gas associated with a galaxy. The major¬ 
ity of the low column density (towards IVhi = 10 13 ° cm -2 ) 
regions reside outside the virial radius of each cluster, and 
primarily towards 2 R v ir- It also appears that the low column 
density HI generally avoids any cluster regions above a few 
million degrees. Although these projections give a nice qual- 
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Virgo-like 


Coma-like 



Figure 1. Density weighted projections of the gas density (top), temperature (middle) and unweighted IVhi (bottom) for the Virgo-like 
and Coma-like galaxy clusters. Note the difference in axes scales between clusters. All projections represent a cubical volume centered 
on the cluster center, 5.0 R v ir on a side. The dashed and dot-dashed circles represent 1 and 2 R v i r respectively. 


itative illustration of the cluster morphological properties, 
they are, of course, not directly observed. We resort, then, 
to synthetic observations of our galaxy clusters in order to 
account for the observability of the HI gas seen in Figure [I] 
and thus see how it would reveal itself in a real observation. 


We discuss the methods of our synthetic observations in the 
following section. 
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2.2 Synthetic Spectra 

As discussed, Y12 used Lya absorption spectra from back¬ 
ground QSO’s to study the warm gas content in the Virgo 
cluster. They compared the inferred neutral hydrogen col¬ 
umn densities and covering fractions to a simulation by cal¬ 
culating the total integrated HI column density along lines 
of sight in a simulated galaxy cluster. In order to make a 
more robust comparison between simulation and observa¬ 
tion, we expand upon this analysis by producing synthetic 
absorption spectra along lines of sight through the simulated 
galaxy clusters. The following section discusses our method 
for generating synthetic spectra which takes into account 
the thermal Doppler broadening of lines as well as the pro¬ 
jected line of sight velocity of the absorbing gas. We then 
automatically fit Voigt profiles with a procedure described 
in Sec. 12331 


to the pure spectra. Although this noise component is simple 
compared to the complicated systematics involved in actual 
COS spectra, it is sufficient to allow a comparison of well- 
resolved column density distributions. We specify a desired 
signal-to-noise ratio (SNR) of 20, slightly better than the 
quality of observations made in Y12. Modifying the SNR 
primarily affects the ability to identify low column density 
features. Changing the SNR does not significantly affect the 
results for absorbers above Ahi ~ 10 13 '° cm -2 . 

Figure [ 2 ] gives an example of the pure spectrum (blue) 
generated using this procedure on the top left, with the noisy 
spectrum (blue) given in the top right. As shown, there are 
two absorption features in the pure spectrum, the smaller of 
which is masked completely by the included noise. Obtaining 
the shown fit (red) is discussed in Section [2. 2. 2| 


2.2.1 Generating Absorption Spectra 


We use the yt toolkit (Turk et al. 2011) to generate absorp¬ 


tion spectra from the simulated galaxy clusters. A complete 
discussion of these methods can be found in [Egan et al.j 
(2014). We outline the procedure below. 


Each absorption spectrum is generated along a line of 
sight through a galaxy cluster at a specified projected dis¬ 
tance from the cluster center (or impact parameter, p). The 
lines of sight generated are all perpendicular to the clus¬ 
ter x-y plane (in the same viewing direction as Figure [I]); 
there is no overlap between lines of sight. We investigated 
the effects of changing viewing angle on the results, and the 
differences are insignificant for the global absorber proper¬ 
ties. However, projection effects do become important in the 
observed velocity distribution of absorbers, as discussed in 
Sec. |4.3| For each cell along the line of sight, the HI column 
density is calculated as Nhi = nm x l, where nm is the HI 
number density and l the length of each cell. The thermal 
Doppler broadening parameter is calculated with the gas 
temperature in the cell, T, as b = \J 2 keT/mH, where kn is 
Boltzman’s constant, and iuh is the mass of the Hydrogen 
atom. Observed absorption of Lya will manifest itself as a 
Voigt profile, having a normalized flux given as F(t) = e ~ T , 
where r is the line optical depth. 

A Voigt profile is computed for each cell along the line 
of sight. Absorption occurs at the natural line wavelength 
of Lya in the frame of reference of the absorbing gas. Thus 
redshift effects and line shifts due to the 3D velocity of the 
gas projected along the line of sight are taken into account. 

For each cluster we observe a cubic volume 5 R v ir on 
a side, centered on the cluster center. This corresponds to 
volumes of 791.5 Mpc 3 and 3375 Mpc 3 for the Virgo-like 
and Coma-like clusters respectively, closely matching the 
observed volume in Y12 and the ongoing Coma observa¬ 
tions, which extend out to p = 2.19 R v n from the cluster 
center. A 50x50 grid of sight lines is placed on each clus¬ 
ter. We generate spectra with a constant wavelength reso¬ 
lution, comparable to the COS G130M grating resolution 
at A=1300 A(used in Y12). The G130M grating operates 
over the wavelength range 1150-1450 A and has a resolv¬ 
ing power (A/AA) R ~ 16,000-21,000. This corresponds to 
a bin width (AA) of 70 mA. To better compare our syn¬ 
thetic observations to COS, we add Gaussian random noise 


2.2.2 Voigt Profile Fitting 

For each spectrum, regions of absorption are identified us¬ 
ing an algorithm adopted fromj Lanzetta, Turnshek fc Wolfe] 
(19871 outlined below. An equivalent width for the absorp¬ 
tion in a given wavelength bin, or e;, is calculated using the 
bin size, AA, flux, Fi, and continuum flux C; as 


e* « AA (1 - Ft/Ci). 


(1) 


Since the generated spectra are normalized to the contin¬ 
uum, the continuum flux, Ci, is unity for each bin. Using the 
uncertainty in each bin, denoted as , a standard deviation 
is computed for each equivalent width as cr ei = AXaF i /Ci. 
Using a passband of 5 bins, a total equivalent width (E;) is 
calculated for a given bin by summing the equivalent width 
of that bin and the 4 adjacent bins. The associated uncer¬ 
tainty is calculated from summing the individual variances. 
Thus, 


i+2 


E >= E 


and 


2 

a Ei 


i+2 

E 


( 2 ) 


A region is flagged as an absorption feature if each bin 
in that region has a total equivalent width larger than the 
associated uncertainty by a specified significance. A spectral 


region is considered a feature if E; > 5a (Lanzetta, Turn 


shek & Wolfe 19871. Each identified absorption feature is 
fit to a Voigt profile with three free parameters: Ahi, b, and 

z. 

Before fitting a given feature, we make initial guesses 
for the absorption region utilizing the procedure discussed 


Dave et al. (1997). Choosing reasonable initial guesses 


is important for obtaining a realistic fit to the absorption 
feature. Each absorption feature is first smoothed using a 
Gaussian filter. A line is placed at the minimum flux of the 
smoothed profile, or the center of a saturated region for satu¬ 
rated lines. For each unsaturated line, Ahi and b are iterated 
over from larger to smaller values until the guessed line is 
above (less absorption) the smoothed feature. For saturated 
lines, Ahi and b are iterated over until the guessed line is 
above the smoothed feature at the edges of the saturated 
region. 

A fit is found from the initial guesses by minimizing x 2 , 
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Figure 2. Top: Example spectra made along a line of sight through the Virgo-like cluster, showing spectra with (right) and without 
(left) noise included with a SNR of 20.0. As shown, the noise obscures the small feature visible at 1218.2 A. Bottom: The JVhi in each 
simulation cell along the line of sight (points), given in velocity space (left), which is the combination of the Hubble flow velocity (vh) 
and projected kinematic velocity, and in the de-projected physical space (right), in terms of just the Hubble flow velocity. The points are 
colored by vh to identify where a cell shifts from the left to right plots. These plots illustrate how the feature fitted in the top spectra 
is matched with a cell in the simulation, as explained in Sec. |2.3| The dashed line and black diamond give the cz and TVhi of the fitted 
absorber in the top right plot. In the bottom plots, vh for each cell is shifted by the cluster Hubble flow velocity such that cells at vh 
= 0 are the same distance from the observer as the distance from observer to cluster center. 


given as 


x 2 = £ 


(Fi-F^y 


(3) 


where F, is the flux in a given bin, is the fitted flux, 
and ap. is the variance given by the uncertainty in each bin. 
This non-linear least squares problem is solved using an im¬ 
plementation of the Levenberg-Marquardt algorithm in the 


SciPy library (Jones et al. 2001). Once a line is fit, the 


resulting fitted flux is subtracted away from the spectrum. 
Then the region is passed through the same absorption fea¬ 
ture identification algorithm discussed above. If there is still 
identifiable flux remaining, a new line is added to the feature 
and the fitting procedure repeats. This is done either until 


no identifiable flux remains, or until the addition of a new 
line worsens the resulting fit. We consider an ideal fit as one 
with a reduced value close to unity. Once lines are fit to 
all of the individually identified features, all of the lines are 
refit simultaneously to the full spectra in a final fitting step 
in order to improve fits on the transitions between adjacent 
lines, especially in very blended features. 

Occasionally, the automatic fitting procedure has trou¬ 
ble separating blended features, attempting to fit the entire 
feature (incorrectly) with a single line. To mitigate this prob¬ 
lem, we have included a step to identify and separate mul¬ 
tiple peaks in an identified absorption region. This method 
is able to properly account for any absorption feature found 
at redshift zero. At high redshift, features become extremely 
blended, with potentially several individual components in a 






































Warm Gas in Galaxy Clusters 7 


given complex. Our method currently has trouble for lines at 
z > 2, requiring manual intervention. However, this redshift 
regime is beyond the scope of this study. 


2.3 Identifying Absorbing Gas in the Simulation 


In addition to making accurate comparisons to observation, 
we are motivated by the ability to use the simulations to 
better understand the physical environment and 3D distri¬ 
bution of the observed warm gas component. We do this 
by associating each identified absorber with a physical re¬ 
gion in our galaxy cluster simulations. Doing this, however, 
has two potential issues: 1) multiple cells along the line of 
sight contribute at various levels to the observed absorption, 
and 2) in reality, the HI absorber seen in spectra depends 
strongly on the geometry and orientation of the HI region. 
This creates an issue in trying to associate a structure/region 
in physical space as an absorber from an identified spectral 


feature (see Egan et al. (2014) for a discussion of this is¬ 


sue). However, at the spatial resolution used here, we find 
that a single cell in the simulation contributes the majority 
of the observed HI column. Although several cells may ul¬ 
timately constitute the physical absorber, the environment 
of these cells does not vary dramatically from the cell with 
most of the HI content, therefore (for simplicity) we only 
associate an identified absorption feature with a single cell 
in our simulations, calling this the physical absorber. 

The absorber(s) identified in a given spectra are located 
along the line of sight in velocity space, rather than physical 
space. The velocity space is the sum of the Hubble flow veloc¬ 
ity of the gas (which is directly related to distance along the 
line of sight from the observer) and the peculiar gas velocity 
projected along the line of sight. Because of this, translation 
to the physical position of the absorbers is impossible with¬ 
out knowing the peculiar gas velocity. This is not possible in 
real observations, but is possible from the simulations since 
we have complete knowledge of the gas motions. 

The process of identifying the gas motions is shown in 
the bottom of Figure[2] On the left, we plot the column den¬ 
sity of each grid cell along the line of sight as a function of 
the line of sight velocity of the gas, given as vi os = vh+Vc, 
where uh = HqD is the Hubble velocity of the gas in a cell 
at distance D along the line-of-sight, and u pec is the gas pe¬ 
culiar velocity projected along the line of sight. This plot 
shows where all of the gas would appear to the observer if 
it was detectable with spectral observations. The spectrum 
on the top right, then, is the sum of the Voigt profiles asso¬ 
ciated with the points on the lower left plot. The Doppler 
shift formula gives a one-to-one correspondence between the 
horizontal axis in the lower left plot, and the wavelength in 
the spectrum in the upper right. The black diamond in the 
lower left plot gives the fitted column density and Doppler 
shifted velocity of the feature shown in the spectrum. As 
this plot illustrates, one cell dominates the total HI column 
density in the region, even though the feature is the sum of 
many cells. 

Since we have full kinematic information in the simula¬ 
tion, we can calculate the component velocities, wh and u pe c, 
separately. We can then de-project the gas by subtracting 
Vpec from i?ios■ This gives the plot on the bottom right, where 
the gas cell column density is plotted as a function of the 
Hubble flow velocity, vh; in other words, this panel simply 


shows the physical location of the gas (expressed in terms of 
the Hubble velocity). Both of these plots are colored by the 
Hubble flow velocity in each cell to readily match the cells 
between the two plots. Comparing these two plots shows just 
how much projection effects shift the gas in velocity space. 
As shown, although the absorber responsible for the feature 
in the spectrum at the top right has a total Doppler shifted 
velocity just below -600 km s -1 (far to the left of most of the 
gas cells in the lower left plot), it is physically located near 
the cluster center (vh ~ 0, lower right plot), and is therefore 
in the middle of the line of sight through the cluster. 


3 PROPERTIES OF WARM GAS IN GALAXY 
CLUSTERS 

We study the warm gas content of two simulated galaxy clus¬ 
ters from the point of view of observable Lya absorbers. We 
discuss the spatial distribution of our cluster warm gas and 
the identified absorbers in Sec. |3.1[ the kinematic properties 
of the absorbers in Sec. |3.2[ and the thermal properties in 
Sec. 13.31 


3.1 Spatial Distribution of Warm Gas 

3.1.1 Cluster Gas Mass Fraction 


We begin by characterizing the diffuse gas around our sim¬ 
ulated clusters. Hot gas (T > 10 G K) is the largest baryonic 
component of galaxy clusters, comprising ~80% or more of 
the cluster baryon mass (VoiT| 2005]). For non-cool core clus¬ 
ters, the average ICM temperature generally decreases ra¬ 


dially outward from cluster center (De Grandi & Molendi 
20021. In the top of Figure[3] we plot the radial profile of the 


gas mass fraction for each of our simulated clusters, sepa¬ 
rated into the three temperature regimes: cold (T < 2 x 10 4 
K), warm (2 x 10 4 K < T < 10® K), and hot (T > 10 6 
K). Here we would like to only examine the diffuse gas as¬ 
sociated with the cluster ICM and surroundings, and ignore 
gas bound in galaxies by placing a cutoff at n < 0.1 cm -3 
for the profiles given in Figure [3] This serves primarily to 
smooth out the cold gas mass fraction curve by reducing 
fluctuations caused by the presence/absence of galaxies in a 
given radial bin. We note that we do not mask this gas in 
our synthetic observations. 

For each cluster, the hot gas dominates (above 80%) 
within the virial radius, and decreases gradually to about 
30 - 40% at the farthest distance from the cluster shown 
(5 R v ir ) ■ Over this range, the warm and cold gas fractions 
increase from, at most, a few percent, to about 30% at 5 
Rvir- The overall trend of the gas fractions is similar in the 
two clusters, but shows a few differences. The Coma-like 
cluster has more warm gas (and less hot gas) around 1-3 
Rvir. It also has a generally lower fraction of cold gas than 
the Virgo-like cluster at all radii. 


3.1.2 Projected map of absorbers 

We now focus our discussion of the cold/warm gas content 
of our simulated galaxy clusters in terms of what is actually 
observable through absorption. We show the projected dis¬ 
tribution of identified Lya absorbers for the Virgo-like (left) 
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Figure 3. The radial profile of the diffuse gas mass fractions for the Virgo-like (left) and Coma-like (right) clusters, as separated into 
three temperature regimes, cold (blue, dash-dotted, T < 2 X 10 4 K), warm (yellow, dashed, 2 X 10 4 K < T < 10 6 K), and hot (red, solid, 
T > 10 6 K). 


and Coma-like (right) galaxy clusters in Figure [4] A circle is 
plotted to represent an identified absorber along a given line 
of sight, with the circle size indicating the column density 
and the color the observed velocity, normalized to the clus¬ 
ter redshift. If a given line of sight has multiple absorbers 
within a column density bin, the circle is subdivided evenly 
among those absorbers. 

One can immediately see the correspondence between 
Figure HI and the total HI column density projections in 
Figure]^ In both clusters, there are several regions of high 
column density seen in the projected map that also appear 
as absorbers in the synthetic observations. In addition, the 
projected map clearly traces filamentary structures at larger 
distances from the cluster center (outside and around 1-2 
Rvir); these are visible also in Figure[l] There is a noticeable 
velocity gradient along groups of absorbers in each cluster. 
The gradients and proximity of absorbers to each other in 
this case immediately suggest a coherent flow or complex 
of absorbers, with a velocity gradient created by a combi¬ 
nation of a true physical velocity gradient and projection 
effects (i.e. coherent flow at an angle to the line of sight). 
An obvious example is the green-yellow structure towards 
the left hand side of the Virgo-like cluster, located mostly 
between 1 and 2 R v i r , but extending out past 2 R v ir- 


3.1.3 Column Density Distribution 


We next turn to a more quantitative analysis of the ab¬ 
sorbers. The column density distribution (CDD) function, 
f(A m ), is defined as the number of absorbers per column 
density and per redshift interval, or f(Am) = d 2 JV/dNdz. 
The column density distribution is often used to charac¬ 
terize Lya absorbers and the associated warm gas in the 
IGM (e.g., |Penton, Stocke fe Shull| |2000[ |Danforth et al.| 
20141. Environmental effects in a galaxy cluster, including 


the interactions with the hot cluster ICM, should result in 
a noticeable departure from a typical IGM distribution of 
Lya absorbers at z = 0. We plot Amf(Am) in Figure [5] in¬ 
stead of f(JVm) to emphasize small differences between the 
Virgo-like cluster (red) and Coma-like cluster (black) distri¬ 
butions. The turnover of the CDD at low column densities is 


an artifact, corresponding to where noise begins to dominate 
over the Lya absorption signal. 

Overall, the Coma-like cluster CDD lies below the 
Virgo-like cluster at almost all column densities. The two are 
comparable only at low column (around Ahi ~ 10 13 cnW 2 ) 


and at high column (around Ahi 


10 1 


highest observed column densities, Ahi > 10 19 cm 2 , both 


cm - ). At the 
cm" 

distributions suffer from low number statistics. Examining 
Figure [l] it is apparent that the Virgo-like cluster contains 
much more absorption at and above Ahi > 10 14 cm" 2 due 
to large area coverage of long, extended HI features. These 
features reside primarily within R v i r of the cluster (in pro¬ 
jection), and correspond to long HI tails that can originate 
from galaxy-galaxy interactions in the cluster, or the ram- 
pressure stripping of gas, which can create tails observable 
in 21cm emission up to 100 kpc ([Oosterloo fe van Gorkom, 


20051. The observable lifetime of these tails in emission (21 


cm and Ha) can be up to at least 500 Myr (Tonnesen & 
Bryan| 20101. Although there has been no study of their 


observability in Lya absorption, presumably the observable 
lifetime of these tails in absorption will be even longer as 
the detectable column densities are substantially lower: Ahi 
~ 10 13 cm" 2 , as opposed to Ahi ~ 10 19 cm -2 . These fea¬ 
tures are absent in the more relaxed Coma-like cluster, per¬ 
haps because these tails have already been heated and mixed 
with the hot ICM. 

We can approximately characterize the CDD as a power 
law with f(Am) oc N^j. We fit this (dashed lines in Figure[5j) 
to each of our distributions over the range 13 < log(Am) < 
14. For the Virgo-like and Coma-like clusters, (3 = -2.284 and 
P = -2.563 respectively. This slope is significantly steeper 
than that observed in the z = 0 IGM, or j3 — -1.68 ± 0.02 
( Danforth et al.| 20141 (we note this slope was fit over the 
range 12.2 < log(Am) < 17). 


3.1.4 Number Density of Absorbers 

The observations for the Virgo cluster and of our simulated 
clusters both show that low column density absorbers gen¬ 
erally avoid the hot X-ray gas and the central R v i r of the 
galaxy cluster in projection. However, we would like to un- 
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Figure 4. Projected map of Lyct absorbers identified in our simulated galaxy clusters, Virgo-like (left) and Coma-like (right). Each 
identified absorber is represented by a circle with size corresponding to the identified column density in one of four IVhi column density 
bins: 10 12-13 , 10 13—15 , 10 15 ^ 17 , and greater than 10 17 cm -3 . We note that gas at HI column densities around 10 12 -10 13 cm -2 is at and 
just below the detection limit in the current observations of the Virgo and Coma clusters. The absorbers are colored by their observed 
velocity, normalized to the cluster redshift. Plotted circles are subdivided if more than one absorber is present in a column density bin. 
In each cluster, the dash and dash-dotted lines denote 1 and 2 R v ; r . These absorber maps have the same orientation as the projection 
in Figure^ 



Figure 5. The column density distribution of identified absorbers 
in the Virgo-like (black) and Coma-like (red) clusters given as 
Nm f(A) 11 ). where f'( JVh [) = 8 2 AT/dNdz. The power law fits over 
the range 13 < log(TVni) < 14 are given for each cluster as the 
dashed lines, with slopes (ft) shown in the legend. 

derstand the distribution of absorbers as a function of the 
physical distance from cluster center. In order to address bi¬ 
ases in our synthetic observations, in Fig [6] we compute the 
three dimensional number density of absorbers in radial bins 
for various absorber column densities. The number density 
is given in units of R~ ; 3 , as opposed to Mpc~ 3 , to account 
for the fact that we scale our sightline grid in terms of R v ir, 
and to readily identify if properties scale with the cluster 
virial radius. We compute the volume probed within a given 


radial bin by summing (xdd 2 for all sightlines, where l is the 
length of a given sightline that falls within a given radial bin, 
and d<5 is the spacing between sightlines (a constant). This 
avoids the problem of taking volume in the i th bin as propor¬ 
tional to R 3 +1 - R 3 , as the number density would then scale 
with the number of sight lines. One could compute the vol¬ 
ume by adding up the volume of “observed” grid cells within 
the radial bins, but this introduces an under-sampling bias 
in the central cluster regions, where the AMR has created 
many more grid cells than in the cluster outskirts. 

Figure [6] shows the number density of all absorbers 
(blue, solid), and absorbers separated into the same column 
density bins as used in Fig. [4] In both clusters, the greatest 
number density of absorbers lies well outside the cluster re¬ 
gion, into the IGM, where the absorber count is dominated 
by low column density absorbers. Outside about 1.5 R v ir for 
each cluster, the number density trends are very similar. By 
this point, the number density of high column absorbers, 
Nhi > 10 15 cm -3 , has decreased (this is more dramatic for 
the Virgo-like cluster), and the number density of low col¬ 
umn absorbers begins to dominate, increasing to the edge of 
the sampled region. High column absorbers show no strong 
trend with radius beyond R v i r for the Coma-like cluster, al¬ 
though there is an uptick beginning around 2.0 R v ir for the 
Virgo-like cluster. This is likely due to the infalling substruc¬ 
ture visible towards the lower right of the Virgo-like cluster 
in Figure [I] Each cluster has surprisingly similar total ab¬ 
sorber number densities outside R v i r although the number 
density is slightly higher in the Coma-like cluster. 

The behavior of low column density absorbers within 
Rvir in each cluster is similar, dropping towards the cluster 
center, although the drop is more dramatic in the Coma-like 
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Figure 6. The number density of absorbers (in units of inverse cubic virial radii) as a function of distance from the cluster center 
for the Virgo-like (left) and Coma-like (right) clusters. The distribution is given for all absorbers (blue, solid), and those with column 
densities of 12.0 < log(Afui) < 13.0 (green, dashed), 13.0 < iog( A) 11 ) < 15.0 (red, dot-dashed), 15.0 < log (A) 11 ) < 17.0 (cyan, dotted), 
and log(Afjji) > 17.0 (purple,solid). 


cluster. In each case, there are no absorbers with column 
densities at A r ni< 10 17 cm -2 found within 0.5 R v ir- In fact, 
for the Coma-like cluster, there are no absorbers at all within 
0.5 Rvir- 

It is the absorbers at Ahi > 10 15 cm -3 that make the 
clusters very distinct. The Virgo-like cluster shows an en¬ 
hancement of these absorbers, such that, while the number 
density for the Coma-like cluster decreases towards cluster 
center, the number density for the Virgo-like cluster actually 
increases within R v ir- This is consistent with the idea that 
higher column density absorbers are associated with galax¬ 
ies, and should be preferentially located within R v i r . Rvir 
also corresponds roughly to where galaxy gas loss effects 


from ram-pressure stripping become significant (Tonnesen, 


Bryan & van Gorkom 2007). This may explain why the 


Virgo-like cluster shows an enhancement in absorbers within 
Rvir, and the Coma-like cluster does not. As the Virgo- 
like cluster is more dynamically active, the ram-pressure 
stripped gas may not have yet had time to thermalize with 
the cluster ICM. This gas is obvious in Figure[l]as long tails 
of HI gas. However, without a greater sample of simulated 
clusters and without better resolution of the ram-pressure 
stripping, it is hard to generalize this apparent difference in 
morphology. 


3.2 Kinematical Properties 

Next, we examine the velocity distribution of the absorbers. 
This is shown in Figure[7]for three different projections. The 
dashed lines in each panel give the dark matter velocity dis¬ 
persion for each cluster as computed directly from the simu¬ 
lation dark matter particles. In most cases, the distribution 
of absorbers is not single-peaked, as would be expected for 
galaxies in galaxy clusters (and for dark matter particles). 
Instead, there is (usually) a clear bi-modality in the velocity 
distributions, with the trough between the peaks centered 
approximately on the cluster systemic velocity, indicating 
that the absorbers are avoiding the central cluster region 
in velocity space. Strictly speaking, however, this does not 
mean that the absorbers avoid the cluster center in phys¬ 


ical space, since the distribution of velocities is subject to 
projection effects. 

To investigate the influence of projection effects, and 
the origin of the troughs in Fig. [TJ we plot the peculiar ve¬ 
locity projected along the line of sight (u pec ) as a function 
of the Hubble flow velocity for each absorber (wh) in Fig- 
ure [8] []] Again, the total observed line of sight velocity of 
any absorber is i>i os = v pec + uh- We center «h on the cluster 
frame, with the vertical line at uh = 0. 

The plots in Figure [5] are subdivided into four regions 
by the vertical black line at vh = 0, and the diagonal black 
line at u pec = -«h- Absorbers with a negative uh are physi¬ 
cally located in front of the cluster, relative to the observer, 
and absorbers with a positive wh are physically located be¬ 
hind the cluster, relative to the observer. Since the observed 
velocity of each absorber is the sum of its value on the verti¬ 
cal and horizontal axes (ui os = vh + u pe c), the diagonal line 
gives a fixed observed velocity of ui os = 0; if Ui os = 0, the 
absorber would appear at cluster center in velocity space. 
This means that absorbers in the upper left quadrant are 
physically located in front of the cluster (uh < 0), but have 
a positive projected kinematic velocity such that v pec > vh, 
causing these absorbers to appear to be behind the galaxy 
cluster in velocity space (or i>i os > 0). Likewise, absorbers 
in the lower right quadrant are physically behind the galaxy 
cluster, but appear to be located in front of the cluster in 
velocity space (ui oa < 0). The absorbers in the remaining 
two quadrants are physically located on the same side of the 
cluster as they appear in velocity space. 

As shown in Figure [8] the majority of absorbers in both 
clusters are located in the upper left or lower right quad¬ 
rants, indicating that most absorbers have large kinematic 


1 We note that the apparent hard cutoff of absorbers at large ab¬ 
solute Hubble velocities is simply because this corresponds to the 
edge of the sampled volume around each cluster. The difference 
in extent between the Virgo-like and Coma-like clusters along the 
horizontal axis is because the sampled volume was a function of 
the virial radius, with the Coma-like cluster being larger by a 
factor of ~ 1.6. 
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Figure 7. Shown is the observed velocity distribution (r>i os = KH+* P ec) for the identified absorbers in the Virgo-like (left) and Coma-like 
(right) simulated galaxy clusters along the z, x, and y viewing axes (from top to bottom). The z axis is the primary viewing direction 
used in this work. The velocities are shifted to the cluster redshift. The dark matter velocity dispersion for each cluster is given by the 
dashed vertical lines. 


velocities pointed towards the cluster center - they are in¬ 
falling. 

The exact fraction of absorbers in these two regions is 
different for these two clusters. This effect seems to be more 
significant in the Coma-like cluster (80% of absorbers), com¬ 
pared to the Virgo-like cluster (68%). The fraction shows no 
obvious trend in the Coma-like cluster with column density, 
although the fraction increases to ~83% for the Virgo-like 
cluster at Nm > 10 15 cm -2 . 


The inset histograms in Figure [8] of the projected pecu¬ 
liar velocities and the Hubble flow velocities reveal a bimodal 
distribution in the projected peculiar velocities of absorbers. 
This is consistent with our earlier observation that the to¬ 
tal velocity distribution of the identified absorbers is also 
bimodal (see Figure [7|. In fact, since the vh distribution 
for both clusters is roughly uniform over the probed vol¬ 
ume, this indicates that the observed bimodal distribution 
is mostly a result of the projected peculiar velocities: or a 
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Figure 8. Plotted is the peculiar velocity of each identified absorber against the Hubble flow velocity for the Virgo-like (left) and 
Coma-like (right) clusters. Note the difference in scales between the vertical and horizontal axes. Absorbers are colored by their column 
density: log(A' r Hi)= 12-13 (black), 13-14 (blue), 14-15 (green), 15-16 (yellow), and 17+ (red). The horizontal axis is scaled such that 
cluster center has a Hubble flow velocity of zero (vertical line). The total observed t>i oa for each absorber is the sum of its peculiar velocity 
and Hubble flow velocity. Absorbers with an observed i>i oa of zero would lie on the diagonal line and would represent absorbers appearing 
at the cluster redshift in projection. 


direct result of the fact that most of the absorbers are in¬ 
falling onto the cluster. Although absorbers do tend to avoid 
the inner cluster ICM (see Sec. 3.1.41, the troughs in Fig¬ 
ure [7] are not indicative of this (since there is no trough in 
the vh histogram of either cluster), but rather are a direct 
result of projected peculiar velocities which change depend¬ 
ing on the orientation of the warm gas filaments relative to 
the observer. It is for this reason that the troughs can shift 
far from vh = 0, and even disappear, in Figure [7] 


These results have all been based on velocities along a 
line-of-sight (i.e. relative to the observer, not the cluster); 
however, we can confirm the infall scenario by looking di¬ 
rectly at the velocity of the absorbers relative to the cluster 
center (v r ). This is shown in Figure [9] It is clear that a 
majority of the absorbers in each cluster represent infall at 
high velocities onto the cluster, as virtually all absorbers 
have large, positive v r . As expected from its larger mass, 
the Coma-like cluster absorber v r distribution generally has 
a higher mean velocity and larger dispersion than the Virgo- 
like cluster. 


Examining the distribution broken down by cluster dis¬ 
tance in Figure |9j there is a noticeable trend: in general, 
moving towards the cluster center shifts the distributions 
towards the right (higher velocities). This is most obviously 
seen going from R v i r > 3 to 2 - 3 R v ir (purple to cyan) in 
both clusters. The shift continues going to 1 - 2 R v ir, but 
less obviously. The bimodal peaks are still present in the 2 
- 3 and R v ir > 3 bins in the Coma-like cluster, but are ab¬ 
sent for absorbers 1-2 R v i r from the cluster center, in favor 
of an almost uniform distribution spanning 500 km s -1 to 
2000 km s -1 . Interestingly, a second peak appears in the 
Virgo-like cluster distribution for absorbers R v ir > 3 from 
the cluster center, which is somewhat responsible for the 


hump in the full distribution. The motion of the absorbers 
towards the cluster center is consistent with the idea that 
most of the observed warm gas represents first infall onto 
the cluster. If any absorbers were to survive passage through 
cluster center, there should be some subset of the distribu¬ 
tion in Figure [9]with negative v r . There are only a handful 
(< 10) absorbers with negative velocities in the Virgo-like 
cluster, and none in the Coma-like cluster. Therefore, if any 
absorbers do survive passage through cluster center, their 
lifetimes must be short. 


3.3 Thermal Properties 


In our final section on the properties of the absorbers, we 
turn to their thermal characteristics. In Figure [To] we show 
the distribution of the absorbers in a temperature-density 
phase diagram (keeping in mind that this is the tempera¬ 
ture and density of the cell with the largest contribution to 
the absorbing gas). To provide context, we also show the 
distribution of all the gas as a gray scale, weighted by mass. 

Looking first at the distribution of all the gas, we see 


that, as expected for such simulations (e.g., Dave et al. 
1 999] ), gas falls largely into three regions: (i) a low-density 
and low-temperature region, similar to the intergalactic 
medium (ICM), (ii) a high-temperature shocked ICM re¬ 
gion, and (iii) cold, dense gas associated with galaxies. As 
expected, the HI absorbers arise primarily from the first 
phase, with most absorbers below about Ahi < 10 14 cm -2 
arising from gas which shows a clear density-temperature 
relation ( Hui fe Gnedin| 19971. For column densities above 
about lO 1 ' cm -2 , the absorbers arise almost entirely from 
dense, cold gas which we expect to be in or near galaxies. 
Interestingly, the population of intermediate absorbers come 
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Figure 9. The distribution of absorber velocities towards the cluster center (v r ) in our Virgo-like (left) and Coma-like (right) simulations, 
with positive velocities indicating motion towards cluster center. The blue (solid) line shows all absorbers, while colored lines indicate 
the distribution split into radial distance bins, with 0-1 R v ; r (green, solid), 1-2 R v ; r (red,dashed), 2-3 R v ; r (cyan, dot-dashed), and 
R v ir> 3 (purple, dotted). 


from gas with a wide range of properties, including some gas 
which is quite warm (T ~ 10 5 K), indicating that it might 
be interesting to look at the width of the absorber profiles 
(e.g., Richter, Fang & Bryan] 2006). 

The metallicity distribution of the absorbers, shown 
in the bottom panels of Figure [TO] indicates that, broadly 
speaking, the IGM-like absorbers have very low metallicities, 
while the absorbers from cold-dense gas are highly enriched. 
The population of absorbers in the transition region between 
warm, diffuse IGM and cold, condensed gas is more mixed, 
with a wide range of metallicities. However, the wide distri¬ 
bution is limited to above about Z = 0.1 Zq for absorbers 
with Nhi > 10 14 cm” 2 . Absorber metallicity is discussed 
further in Sec. 14.21 


4 DIRECT COMPARISON TO OBSERVATION 

The previous section demonstrates that the majority of ab¬ 
sorbers identified in our simulations come from gas which is 
of low density falling into the cluster for the first time. We 
now turn to making direct comparisons and predictions for 
observations of absorbers in and around galaxy clusters. It 
is important to keep in mind that we currently only have 
the Virgo Cluster observations to compare to directly and 
that the number of sightlines through a cluster will always 
be very limited compared to our simulation results. We be¬ 
gin with the covering fraction, as this can be computed for 
small sample sizes. We then go on to discuss the metallici¬ 
ties and the position-velocity distributions of the absorbers 
in the two simulated clusters. 


4.1 Covering Fraction and Column Density 
Distribution 

The covering fraction of absorbers, f COV er, can be computed 
as the fraction of sight lines that contain a feature above a 
given column density. We present the covering fraction for 
our Virgo-like and Coma-like clusters in Figure EU along 



log Nhi cm 2 


Figure 11. Comparison of covering fraction as a function of Nhi- 
The covering fraction is defined as the fraction of lines of sight 
with identified absorption at or above the given Nhi- The syn¬ 
thetic observations are shown with solid lines for the Virgo-like 
(black) and Coma-like (red) clusters, with the total (integrated 
through all cells along the line of sight) Nhi for each shown 
with dashed lines. These are compared to the Virgo observa¬ 
tions (black, dash-dot). The shaded region marks the complete¬ 
ness limit at log( Nhi cm -2 ) < 13.3 given in Y12. 


with a direct comparison to the actual Virgo covering frac¬ 
tion from Y12. The gray region in this figure denotes the 
completeness limit of Vhi> 10 13 ' 3 cm -2 given in Y12. The 
covering fractions for our Virgo-like and Coma-like galaxy 
clusters are similar with each other. Yet, above the complete¬ 
ness limit the Coma-like cluster has a lower f CO ver than the 
Virgo-like cluster, in agreement with the column density dis¬ 
tribution in Figure [5] Both simulated clusters, however, lie 
well below the actual covering fraction of the Virgo cluster 
at all column densities. The discrepancy varies by a factor 
of up to 7 at Nm < 10 15 cm -2 , to about 3 at JVhi > 10 15 
cm” 2 . 

To ensure that our synthetic observation methods are 
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Figure 10. Phase plots of all the gas (greyscale) and absorbers (colored circles) for the Virgo-like cluster (left) and Como-like simulations. 
In the top set of panels, the absorbers are color coded by column density, while in the bottom, they are color-coded by metallicity. 


not the source of this disagreement, we plot the covering 
fraction of the total HI content in the simulation integrated 
along each line of sight as the dashed lines in Figure im 
The integrated HI covering fraction reproduces the relation¬ 
ship between the two clusters in that the Virgo-like cluster 
generally has more absorption than the Coma-like cluster 
(above the completeness limit). In addition, it reproduces 
the larger covering fraction of the Coma-like cluster at low 
column densities (below the completeness limit). Above the 
completeness limit of Y12 (which we tried to reproduce us¬ 
ing noise comparable to the COS observations of Y12), the 
integrated f COV er corresponds well with the synthetic obser¬ 
vations. Therefore, the missing HI content is not an artifact 
of the synthetic observational technique, but corresponds to 
a real lack of HI in the cluster simulations. 

To examine morphological differences between the two 
simulated clusters, and further compare to observations, we 
tabulate the covering fraction for each cluster at various col¬ 
umn densities separated into R v ir bins in Table [2] The col¬ 
umn densities correspond to those chosen in the Y12 Virgo 
results, corresponding to equivalent widths of W = 65, 100, 
150, 200, and 300 mA, assuming a Doppler broadening of 
b = 30 km s -1 . The disagreement between the Virgo cluster 
and our synthetic clusters is clear here, with covering frac¬ 
tions of unity at Nm= 10 13123 cm~ 2 in Virgo, and fractions 
ranging from several to 20% in our simulated clusters. The 


general trend seen in the Virgo observations is that covering 
fraction decreases with increasing column density and in¬ 
creases away from the cluster center. This trend is seen also 
in our Coma-like cluster, and extends outside 2 R v ir, a region 
not probed by the Virgo observations. However, this behav¬ 
ior is not matched in our Virgo-like cluster at A1 hi> 10 13 ' 572 
cm~ 2 , as the maximum covering fraction at these column 
densities occurs within the virial radius of the cluster (in 
projection). At all column densities, however, the fraction 
at 2 - 2.5 Rvir is greater than that at 1 - 2 R v ir- 

Consistent with what is shown by the covering fraction, 
the column density distribution of our simulated clusters lies 
below that given in Y12, as shown in Figure[l2] We plot here 
the normalized probability density function, rather than the 
CDD given in Figure[5] The Virgo-like cluster contains more 
absorbers compared to the Coma-like cluster except at very 
high and low column densities. The largest disagreement 
between the simulated clusters and Virgo is, again, near 
Nm= 10 14 cm~ 2 , though the disagreement extends over all 
column densities. The limited sight lines through the Virgo 
cluster and small number statistics make rigorous compari¬ 
son to our simulations at A r m> 10 15 cm -2 difficult. 

Finally, the typical number of absorbers per sight line 
in our cluster simulations also lies below that seen in Y12, 
even though in both cases the sightlines probe similar re¬ 
gions (in velocity space). In the Virgo cluster observations 
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Virgo 


Virgo-like 



Coma-like 


JVhi 

0-1 R v i r 

1-2 R vi r 

0-1 Rvir 

1-2 Rvir 

2-2.5 Rvir 

0-1 Rvir 

1-2 Rvir 

2-2.5 Rvir 

13.123 

1.00 — 0.15 

1.00-0.17 

0.086 

0.134 

0.159 

0.065 

0.175 

0.202 

13.338 

0-60i°i 3 

1.00-0.14 

0.083 

0.093 

0.106 

0.033 

0.090 

0.102 

13.572 

0-50t°; 4 ® 

0.86±Hs 

0.073 

0.059 

0.072 

0.027 

0.052 

0.050 

13.758 

0.27±°i® 

0.75°j 41 17 

0.074 

0.037 

0.051 

0.024 

0.035 

0.034 

14.109 

o-osi °;£ 2 

o-se ±8;1| 

0.067 

0.030 

0.033 

0.017 

0.018 

0.022 


Table 2. Covering Fractions for the Virgo-like and Coma-like clusters in selected IVhi and R v j r ranges. 



Figure 12. The normalized column density distribution for the 
Virgo-like (black) and Coma-like (red) clusters in comparison to 
the observed distribution in the Virgo cluster from Y12 (black, 
dash-dotted). The gray region marks the completeness limit of 
Y12 given as 1 Vhi = 10 133 cm -2 . 


from Y12, 70% of COS sightlines (6 of 9) had more than 
one absorber, and two sight lines contained five absorbers 
each. The Y12 STIS and GHRS sightlines also show this 
percentage of multiple absorbers. From the simulations, we 
find only 22.6% of the sightlines with absorbers have more 
than one in the Virgo-like cluster and 24.7% in the Coma¬ 
like cluster. In addition, the maximum number of absorbers 
per sightline is only 3 (4.7% in the Coma-like cluster, none 
in the Virgo-like). The differences between the Virgo obser¬ 
vations and our simulated clusters show that the difference 
in absorber frequency is noticeable on an individual sightline 
basis. The differences do not manifest themselves solely in 
the global absorber properties (i.e. Figures 
Tabic [2}. 
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and 
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In our comparison to observations from Y12, we note 
that it is possible Y12 included absorbers not directly associ¬ 
ated with the Virgo cluster or its substructures. We compare 
against the full velocity range of Y12, or 700 < cz < 3000 
km s” 1 ; for reference, M87 is located at about 1307 km s” 1 , 


km s” . As shown in Figure |7| the absorbers in 
lated clusters mostly lie within about ov of cluster systemic 
velocity. Making a 2cr v and 1 <r v velocity cut on the Y12 ab¬ 
sorbers reduces the covering fraction from about 80 % at the 
completeness limit (A r m= 10 13 ' 3 cm” 2 ; see Figure 111, to 
72% and 36% respectively. At Nm= 10 14 cm” 2 , the cover¬ 
ing fraction stays at 32% for the 2ov cut, but goes down to 


and 10% at la v . These velocity cuts would bring the Y12 
observations into much closer agreement to our Virgo-like 
cluster at Nm> 10 14 cm” 2 , but a disagreement by a factor 
of a few remains below Nm= 10 14 cm” 2 . One should keep in 
mind, however, that the Y12 observations are missing part 
of the velocity range of the Virgo cluster. 

We also approach the discrepancy from the simulation 
side and investigate the possibility that the disagreement 
seen in Figures and [T2] is due to the difference in ob¬ 
served and simulated cluster volume. Using H 0 = 100 km 
s” 1 , and ignoring peculiar velocities, the full Y12 velocity 
range corresponds to probing 23 Mpc of physical distance 
along the line of sight in Virgo. Our synthetic observations 
probe 5 R v i r along a line of sight, corresponding to 9.25 and 
15.0 Mpc for the Virgo-like and Coma-like clusters respec¬ 
tively. We repeat our analysis using the same lines of sight in 
each cluster, but increasing their length to probe the entire 
high resolution region in each cluster, corresponding to 45 
and 60 Mpc in the Virgo-like and Coma-like clusters respec¬ 
tively. As expected, the covering fraction and CDD increases 
with the larger observed volume, but the increase lies only 
for absorbers at Ahi< 10 13 cm” 2 , by a factor of about 2.5 
(fcover = 0.33 for the Virgo-like cluster at N m = 10 12 cm” 2 ). 
Increasing the observed volume in the simulation (thereby 
including fore/background structures) does not appreciably 
bring the simulations closer into agreement with Y12. Ulti¬ 
mately, removing fore/background absorbers in observations 
presents a challenge in comparing between simulation and 
observation. Further observations will help clarify if the dis¬ 
agreement seen here is universal, or if it is due to peculiari¬ 
ties of the Virgo cluster. 


4.2 Absorber Metallicity 

Since complete distance and kinematic information of the 
absorbing gas in unavailable in observations, absorber metal¬ 
licity could be a valuable discriminator between the IGM 
and galactic origins of absorbers. Within the spectral range 
of the COS observations in Y12 and the ongoing Coma ob¬ 
servations, the Si III [1206], Si II [1260], OI [1302], and CII 
[1334] lines can potentially be used to obtain an estimate 
of the absorber metallicity. The IGM at z = 0 has a low 


= 544 metallicity, Z < 0.1 Zq (Danforth et al. 

2014 

1, making it 

sirnu- distinct from the typical metallicity of ga 

act ic out Hows, Z 


> 0.1 Zq up to Z ~ Zq (Creasey, Theuns & Bower 20151. 

In our simulated galaxy clusters, we find an important 
distinction in absorber metallicity with both absorber col¬ 
umn density and distance from the cluster center. In both 
clusters, almost all absorbers with column densities above 
about Ahi = 10 14 cm” 2 have Z > 0.1 Zq, and up to a 

































16 A. Emerick et. al. 


few times solar, regardless of their location in the cluster 
environment. High column density absorbers, then, are gas 
processed through galaxies (galactic waste). Absorbers with 
column densities below 10 14 cm -2 span the entire range of 
observed metallicities, from virtually pristine up to solar, 
yet most lie below 0.1 Zq with the greatest concentration 
at virtually pristine metallicities. Low column density gas 
is then easily associated with the IGM, yet some may still 
have been processed in galaxies. 

Metallicity may also prove valuable in disentangling 
projection effects. In projection, absorbers with metallici¬ 
ties between galactic waste and pristine gas are only found 
outside a distance of about 1.5 R v i r - Within this radius, ab¬ 
sorbers appear to either be enriched or pristine (with very 
little middle ground). Metallicity as a function of physical 
distance from the center of the simulated cluster provides 
a stark contrast. We find that all absorbers within a phys¬ 
ical distance of about 1.5 R v i r and 1.3 R v i r for the Virgo- 
like and Coma-like clusters respectively, have metallicities 
greater than Z ~ 0.05 Zq. It is only outside these radii that 
absorbers have low metallicity or are pristine. Therefore the 
infalling IGM absorbers do not penetrate inside the clus¬ 
ter virial radius without first being heated and/or ionized. 
Since the vast majority of the Virgo absorbers are low col¬ 
umn density and no metal lines were detected, the majority 
are likely to lie beyond the virial radius as suggested by Y12. 


4.3 Velocity Distribution of Absorbers 

The velocity distribution of the simulated absorbers in Fig¬ 
ure 0 shows that most absorbers sit in front of or behind 
the cluster in velocity space. Figures [8] and [9] show that the 
majority of the absorbers that appear in front of and be¬ 
hind the galaxy cluster in Figure [7] are in fact physically 
located on the opposite side of the cluster as they represent 
gas falling in along cosmic filaments. Unfortunately, we are 
unable to directly compare our velocity distributions to the 
Virgo cluster observations because the damping wings of the 
Milky Way render a large region of velocity space unobserv¬ 
able in the Virgo cluster (Y12). The observed distribution is 
truncated right around where a forward peak would lie (if it 
exists in the Virgo cluster) making it difficult to interpret the 
nature of the distribution. However, there is a comparative 
lack of absorbers at the Virgo cluster systemic velocity. The 
intuitive explanation for this trough is that, since a majority 
of absorbers avoid the hot ICM, there is a decrement in the 
central cluster regions in velocity space. However, since the 
redshift projection effects of these absorbers are dramatic, 
we examined the observed velocity distribution along two 
other, orthogonal viewing angles (‘x’ and ‘y’) from our pri¬ 
mary observations (along the ‘z’ axis). These distributions 
are included in the bottom two rows of Figure [T] For both 
clusters, all but one projection shows a bimodal distribution 
of absorber velocities. Given the bimodality is not univer¬ 
sal, we conclude that the observed velocity distribution of 
absorbers in a galaxy cluster is biased by the unique na¬ 
ture of the infalling material, namely the location and size 
of infalling filamentary material relative to the observer. 


4.4 Projected Distribution of Absorbers 


Observations can examine the distribution of column den¬ 
sity with projected distance from the cluster center. Y12 
found tentative evidence for a positive correlation between 
absorber column density and projected distance from the 
cluster virial radius when absorbers associated with Virgo 
substructures were excluded. Our Virgo-like cluster does not 
show this behavior in projection (black, dashed) as shown 
in Figure [13] which gives the median absorber column den¬ 
sity in radial bins; although we make no cuts to remove 
absorbers associated with substructures. The median col¬ 
umn is greatest near the cluster center, at Ahi > 10 15 
cm -2 , and decreases towards the cluster outskirts, remain¬ 
ing roughly constant at around 10 13 cm -2 at projected dis¬ 
tances greater than 1.5 R v ir- The radial profile of the Virgo- 
like cluster (black, solid) shows the transition from cluster to 
outskirts more clearly, however, remaining near 10 15 cm -2 
up to R = 1.25 Rvir, then dropping off quickly to Ahi ~ 10 13 
cm -2 at R > 1.5 R v ir- 

The Coma-like cluster radial profile (red, solid) is qual¬ 
itatively similar to the Virgo-like cluster, yet has a signifi¬ 
cantly higher median column within R v i r dropping quickly 
to 10 13 cm~ 2 outside the virial radius. This is a dramatically 
sharper transition than is shown in the Virgo-like cluster. 
Interestingly, the projected profile of the Coma-like cluster 
is flat, at Ahi ~ 10 13 cm -2 , at all impact parameters. In 
projection, the few high column central absorbers in the 
Coma-like cluster are drowned out by the low column fore¬ 
ground/background absorbers. This may be similar to what 
is observed for Virgo, as the distribution of column densities 
with impact parameter is relatively flat before the removal 
of sub-structure absorbers (Y12). 

The Coma-like cluster stands in contrast to the ob¬ 
served relationship of Ahi with impact parameter for galax¬ 


ies (e.g., Prochaska et al. (20111; Stocke et al. (20131; 


Tum- 


linson et al. ( 20131). Ahi has been consistently found to in¬ 
crease with decreasing impact parameter for observed galax¬ 
ies. Although our Virgo-like cluster shows this relationship, 
the flat profile of the Coma-like cluster indicates there may 
be a change in this observed relationship for the most mas¬ 
sive dark matter halos. This suggests there is some definite 
mass scale at which point the nature of gas accretion onto 
dark matter halos changes, and a majority of the accreted 
gas remains hot and unable to cool. This will be interesting 
to test further with future observations and simulations of 
galaxy clusters. 


5 DISCUSSION 

We discuss the implications of our study on the understand¬ 
ing of gas flows in and around galaxy clusters in Sec. |5.1| 
The effects of resolution on studying the Ly a absorbers in 
our clusters is discussed in Sec. |5.2[ Finally, limitations of 
our simulations in capturing physical processes relevant to 
the formation and destruction of absorbers, including cluster 
specific processes, are summarized in Sec. |5.3[ 
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Figure 13. The median column density of absorbers as a func¬ 
tion of radial (solid) and projected (dashed) distance from cluster 
center for the Virgo-like (black) and Coma-like (red) clusters. In 
each case, there is a negative trend with distance. 


5.1 Gas Flows in Galaxy Clusters 


We find that warm, absorbing gas in the environment sur¬ 
rounding galaxy clusters is generally infalling in our simula¬ 
tions. The warm gas identified in absorption that is located 
within the virial radius of the cluster is most likely galactic 
in origin, generally having both a higher observed HI column 
density and a higher metallicity than an absorber located in 
the pristine, diffuse IGM. Outside of the virial radius, the 
total population is dominated by lower column density and 
metallicity absorbers. These are associated with the IGM 
filaments that are feeding the growth of our two simulated 
galaxy clusters. Like the absorbers within the virial radius of 
the cluster, these filamentary absorbers are fast moving, in¬ 
falling material. However, the absorbers at large radii more 
clearly represent coherent gas flows from the IGM. 

Although we only examine two simulations in this study, 
we find that morphology and cluster dynamics are likely to 
affect the resulting presence and distribution of warm, neu¬ 
tral gas. Our more relaxed and massive Coma-like cluster 
exhibits noticeably less absorption within R v i r than the still 
merging Virgo-like cluster. As the Virgo-like cluster relaxes, 
it is likely that many of the galactic absorbers within the 
virial radius will reach thermal equilibrium with the hot 
IGM, relaxing to a state similar to that of the Coma-like 
cluster. Due to its large mass, however, the Coma-like clus¬ 
ter has significantly more filamentary absorption than the 
Virgo-like cluster, cand these filaments have more physical 
volume and contained gas mass. Interestingly, the flat Ahi 
profile for the Coma-like cluster in Fig. [13] suggests that 
there is a dark matter halo mass regime where the nature 
of warm, neutral infalling gas changes, when compared to 
the Virgo-like cluster, and similar observed profiles for Ahi 
around galaxies (e.g., |Prochaska et al.| ( |2011| ) ; |Stocke et al.| 
(20131; Tumlinson et al. (20131). 


5.2 Resolution Study 

The resolution in any simulation is important to consider in 
order to assess the validity of the results. In our case, too 


low a resolution could prevent proper accounting of high 
column density absorbers located in very overdense regions, 
and can spread diffuse absorbers too thin, reducing their ob¬ 
served column density (where again, Nhi ~ (nm where l is 
a given cell length, and nm is the neutral hydrogen num¬ 
ber density). Lukic et al. (20141 found a resolution of 20 
h -1 kpc was necessary for 1% convergence in the Lya forest 
at z = 2. Over the region around each cluster probed by our 
observations, we have a minimum and maximum cell res¬ 
olution of 7.3 and 233 kpc respectively. Averaging over all 
cells, the Virgo-like and Coma-like clusters have a resolution 
of 31 and 34 kpc, respectively. However, a simple extrapo¬ 
lation like this is challenging as the gas over density giving 
rise to a given Ahi absorber is higher at low redshift, due to 


the evolution in the ionizing background (e.g., Dave et al. 


2010 ). 


It is computationally difficult to improve the resolution 
of our cluster simulations; however, we can test the resolu¬ 
tion dependence of our results by examining the convergence 
of a smaller, uniform region. In particular, we examine the 
column density distribution function f(N) in 3 unigrid sim¬ 
ulations of a 20.0 h -1 Mpc box using the same parameters 
and included physics as the cluster simulations described in 
Section[2] Essentially, these simulations probe the z = 0 Lya 
forest in the IGM (sacrificing the higher column density ab¬ 
sorbers which are due to dense gas which is not resolved in 
these unigrid calculations). Although the grid in our cluster 
simulations is certainly more complex, with adaptive mesh 
refinement, these unigrid simulations, performed with 128 3 , 
256 3 , and 512 3 grid cells and dark matter particles, provide 
a straightforward relationship between f(N) convergence and 
grid size in a computationally tractable manner. To test the 
effect of AMR, we conduct a fourth AMR enabled simula¬ 
tion with an initial grid of 128 3 . These simulations have a 
corresponding resolution at z = 0 of 232.9 kpc, 116.4 kpc, 
and 58.2 kpc, with the 128 3 AMR simulation reaching the 
same minimum and maximum resolution as our cluster runs 
(7.28 and 232.9 kpc), with an average resolution of 86 kpc. 

In Figure [l4] we plot Nmf(Nm) for these simulations. 
As shown, the column density distribution is well converged 
for low column density absorbers (Nm < 10 14 cnC 2 ) even at 
128 3 (black) resolution (232.9 kpc grid size). However, the 
distribution is not well resolved at Ahi > 10 14 cm -2 at this 
resolution, and is only converged up to Ahi ~ 10 14 ' 5 at 512 3 
(blue). At column densities above this point, the fixed grid 
resolution of 58.2 kpc is insufficient to completely capture 
these absorbers, and the CDD begins to deviate from a single 
power law. Fortunately, in our cluster simulations, we are 
not restricted to a fixed, uniform grid. 

For the Virgo-like and Coma-like clusters, 77.8% and 
67.2% of the probed volume has resolution equal to or bet¬ 
ter than the 256 3 simulation. Likewise, 31.6% and 21.7% 
of each cluster has resolution at or better than the 512 3 
simulation. Only overdense regions refined by the AMR in 
the cluster simulations reach resolution better than the 512 3 
run, or 5.8% and 3.9% of the total cluster volume. Our clus¬ 
ter simulations are thus well converged for absorbers at Am 
< 10 14 0 cm -2 , and are fairly well converged at Am < 10 14 '' 
end 2 . However, the results (like the 512 3 run) are probably 
not well converged at higher column densities of Am > 10 15 
cm~ 2 . Most of the observed Virgo absorbers are below this 
column density and so overall comparisons between the sim- 
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Figure 14. A resolution study of the column density distribu¬ 
tion given as Njji f(Nn I) for our three unigrid simulations of a 
20 h -1 Mpc box at 128 3 (black), 256 3 (red), and 512 3 (blue). In¬ 
cluded also is the 128 3 run with one level of AMR enabled (green). 
Observational results from |Danforth et al.] ( |2014[ ) are shown in 
purple. 


ulations and observations will not change with resolution. 
However, the specific column densities of the high column, 
galactic absorbers may vary somewhat with increased reso¬ 
lution runs. 


5.3 Simulation Uncertainties 

Although resolution is likely the largest source of uncertainty 
for high column density absorbers, it is not the dominant 
source of uncertainty at modest column densities as shown 
by the convergence in Sec. |5.2| We discuss additional uncer¬ 
tainties and missing physics in our simulations that could 
play a role in creating or destroying warm neutral gas that 
would be visible as Lya absorption. 


5.3.1 Ionizing Background 


The adopted ionizing background plays a significant role in 
establishing the neutral fraction at all redshifts. The ioniza¬ 
tion rate, T, is well constrained at high redshifts (z > 2), 
with excellent agreement between the observed and simu¬ 
lated Lya forest (e.g., Croft et al. 2002 Viel et al. 2005| 


Seljak et al. 2005). However, the low redshift Lya forest 


is poorly constrained, and the z = 0 forest especially, is 
not well understood. Although pre-COS observations have 


placed constraints on the low redshift IGM ( Dave et al. 
[2010, and references therein), COS has only recently allowed 


for the largest sample of low redshift Lya observations (Dan 


forth et al. 2014). There are many simulations of the high 


redshift Lya forest, yet only a few works examine the for¬ 


est in the local Universe (Paschos et al. 2009 Dave et al. 


2010 Egan et al. 2014 Kollmeier et al. 20141. Most of these 
employ the |Haardt fe Madau| ( |2001| ) ionizing background 
and suggest fair agreement with observation, although gen¬ 
erally underreporting the HI column density distribution by 
a factor of a few. To date, only one work has examined the 
low redshift Lya forest using the updated |Haardt fc Madau| 


(20121 background, Kollmeier et al. (20141, which strongly 
suggests that the updated background results in too much 
HI in the IGM at 2 = 0 at low to moderate column densi¬ 
ties. This suggests a significant uncertainty in the sources of 
ionizing background radiation at low redshift, which could 
substantially change the neutral HI content in our simula¬ 
tions. 

As shown in Figure |14[ our highest resoluti on unig rid 
simulation, 512 3 (blue), agrees qualitatively with 


Kollmeier 


et al. (20l4l, with an overproduction of HI compared to 


the observations in Danforth et al. (20141 (purple), at least 
at low column densities. At IVhi > 10 5 cm -2 , our IGM 

simulation starts to drop below the observations, where res¬ 
olution begins to have an effect. 

Although we only ran our cluster simulations with 
the HM12 background, we can analytically estimate how 
a change in the background would affect the observed Nm- 
Since the gas studied here has a large ionization fraction, we 
can assume ionization equilibrium. Considering only elec¬ 
trons and Hydrogen, the Hydrogen neutral fraction, and 
thus the HI column densities, will scale inversely with the 
choice of background ionization rate (T) (see Rauch (19981 
for a complete summary of this argument). For HM01, Thi 
is roughly a factor of 5 higher at z = 0 than in HM12. Using 
the HM01 background instead of HM12, then, would reduce 
the neutral fraction by a factor of 5, effectively shifting the 
lines in Figure [l4| to the left. 


5.3.2 Galaxy Outflows 


It is currently computationally prohibitive to run cosmolog¬ 
ical simulations of galaxy clusters to high enough resolution 
to completely capture the galaxy physics responsible for sub¬ 
stantial outflows (indeed, this is challenging even for zoom 
calculations of individual galaxies). Relevant physics can in¬ 
clude a multi-phase ISM, detailed effects of ram-pressure 
stripping of the ISM, star formation feedback, and AGN 


(e.g., Springel & Hernquist 2003 Bower et al. 2006 Cover 


nato et al. 2007). However, with a maximum resolution of 


several kpc, the galactic outflow processes must be addressed 
with subgrid physics models. As mentioned in Sec. ED our 
simulations include star formation, which contains the “dis¬ 
tributed feedback” model for star formation and supernova 


feedback discussed in Smith et al. (2011). Although this 


model still suffers from overcooling and hence overproduc¬ 
tion of stars and metals, in general [Smith et al. ([2011) shows 
reasonable agreement with observations of the distribution 
of metals in the IGM. 

Although the influence of various feedback models on 
the Lya forest in Enzo simulations has not been tested, 


Dave et al. (2010) discusses the effects of various galac¬ 


tic wind models on their SPH simulations of the low red¬ 
shift IGM. Their preferred wind model, a momentum driven 
wind, produces a distribution of low redshift Lya absorbers 
consistent with the then-available observations. In general, 
they found that turning on their wind model had no signifi¬ 
cant effect on low column density absorbers, those with Nm 
< 10 14 cnU 3 . The differences between wind models became 
significant for higher column density absorbers, especially 


for lines of sight with low galactic impact parameters (Ford 
[et al.[ |2013[ ). It is difficult to judge how much the cluster en¬ 
vironment would modify these results obtained for the IGM. 
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If much of the observed HI absorption in the Virgo cluster is 
galactic in origin, it is likely that feedback processes and the 
effects of the cluster environment on galaxies must be exam¬ 
ined further as possible additional sources of warm HI. As 
demonstrated in Dave et al. (20101 and Ford et al. (20131, 
however, this would likely only serve to increase the number 
of absorbers observed at Ahi > 10 14 cm -2 , and therefore 
would not be able to fully explain the differences between 
our simulations and the Virgo observations. 


5.3.3 AGN and Cool-core Clusters 


Finally, we note that the Virgo cluster is a cool-core cluster 
that exhibits a drop in ICM temperature towards cluster 
center, and strong signs of active AGN feedback. The Coma 
cluster does not have significant current AGN jet heating 
and is not a cool core cluster. The effect of cool-cores on 
the presence of warm gas within the ICM of the cluster is 
unknown, yet cluster properties may be dependent on the 


presence/absence of a cool-core (e.g., Frank et al. 20131. 
The physics required to generate a cool-core cluster in our 
simulations is not included, and is currently computationally 
prohibitive to employ in detail in cosmological simulations, 
although the effects of AGN in clusters has been studied 
before with Enzo ( |Li fc Bryan| [2014|>. Inclusion of AGN on 
cool and warm diffuse gas is likely only to be significant (if at 
all) in the very central regions of the cluster, well within R v ir- 
Our ongoing observations of the Coma cluster will be one 
way to address observationally the comparative differences 
between warm gas in cool core and non-cool core clusters. 


6 SUMMARY 

The flow of gas into and throughout galaxy clusters is fun¬ 
damental to their evolution. Warm gas traces gas flow across 
cluster scales, between galaxies and the ICM, and between 
the cluster environment and the IGM. Our work is the 
first systematic study of the warm gas content in simu¬ 
lated galaxy clusters in the context of Lya absorption. These 
simulations provide a new understanding of the spatial and 
kinematic properties of the absorbing gas, as well as insight 
into their environment. Utilizing synthetic observations to 
carry out this study allows us to make direct comparisons 
to previous observations of the Virgo cluster, and to estab¬ 
lish what we expect from future cluster observations. We 
summarize the main results of this work below. 

(i) In each of our clusters, the warm and cold gas mass 
fraction is small within R v i r (less than a few percent) but 
increases substantially (to about 30% at R > R v ir) transi¬ 
tioning away from the hot cluster environment into the IGM. 
The population of absorbing material traces this gas, gen¬ 
erally avoiding the hot cluster center (with the exception 
of some high column density absorbers, which correspond 
to galactic waste). In particular, we show that the numer¬ 
ous low column density absorbers trace the filaments of the 
cosmic web feeding gas into the cluster. 

(ii) Low column density absorbers dominate the total ab¬ 
sorber counts in each cluster, and dominate the number den¬ 
sity of absorbers in each cluster outside 1 R v ir and 2 R v i r 
for the Coma-like and Virgo-like clusters respectively. The 


increased number density (in units of R v i r -3 ) of low column 
absorbers in the Coma-like cluster indicates greater filamen¬ 
tary feeding in the more massive, more relaxed of our two 
clusters. 

(iii) The absorbers in each of our simulated galaxy clus¬ 
ters are generally fast moving and infalling for the first time 
into the cluster, with peculiar velocities sufficiently larger 
than their Hubble velocity such that they are redshifted to 
the opposite side of the cluster central velocity. This is con¬ 
sistent with the idea that absorbers trace infalling filaments 
which do not survive first infall. 

(iv) In general, the observed velocity distribution of ab¬ 
sorbers is bi-modal, with two peaks separated by a trough 
centered roughly on the cluster systemic velocity. The rel¬ 
ative size of these peaks, and the width of the trough, is 
determined by projection effects. 

(v) When compared to the Virgo cluster observations of 
Y12, both of our simulated clusters show less than the ob¬ 
served warm HI content. In terms of covering fraction, the 
disagreement is a factor of several, with f CO ver ~ 0.4 for Virgo 
and f CO ver ~ 0.05 for our simulations at Nm = 10 14 cm -2 . It 
remains to be seen whether the Lya absorbing warm gas in 
the Virgo cluster is representative of galaxy clusters in the 
nearby universe. Future observations will be important for 
further interpreting this disagreement. 

(vi) Metallicity is an important discriminator of absorber 
origin. High column absorbers have metallicities character¬ 
istic of gas processed in galaxies, and low column absorbers 
are predominately pristine. No low metallicity absorbers are 
found within the cluster virial radius, suggesting that the in¬ 
falling IGM absorbers heat up or ionize quickly, even before 
entering the cluster environment. 

(vii) A resolution study indicates that our low-column 
density population Nm < 10 14 ' 5 cm' 2 is reasonably well- 
converged, but that we probably systematically under pre¬ 
dict higher column-density absorbers. However, this is not 
enough to explain the above difference in f CO ver between 
Virgo and our simulations. 

This work provides important insight into cluster warm 
gas content. Cluster Lya absorbers correspond to fast mov¬ 
ing, filamentary infall, which is enhanced in more massive 
clusters. The presence/absence of absorbers within R v i r is 
dependent upon cluster morphology, and we expect dynam¬ 
ically young clusters to have enhanced absorption within 
Rvir- These simulations can be used to examine the corre¬ 
spondence between galactic absorbers and the cluster galax¬ 
ies from which they came. In addition, higher resolution ver¬ 
sions of these simulations can be used in detail to trace the 
origin and fate of the identified absorbers, examining the 
lifetimes of absorbers with different origins (filamentary vs. 
galactic). Future cluster observations, including the ongoing 
COS Coma observations, will be important to understand¬ 
ing how absorbers vary with cluster mass and morphology 
over time. 
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